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Abstract 

Phenomena analogous to ground state quantum phase transitions have recently 
been noted to occur among states throughout the excitation spectra of certain 
many-body models. These excited state phase transitions are manifested as simul- 
taneous singularities in the eigenvalue spectrum (including the gap or level density), 
order parameters, and wave function properties. In this article, the characteristics 
of excited state quantum phase transitions are investigated. The finite-size scaling 
behavior is determined at the mean field level. It is found that excited state quan- 
tum phase transitions are universal to two-level bosonic and fermionic models with 
pairing interactions. 

PACS: 03.65.Fd, 03.65.Sq, 64.60.-i 



1 Introduction 

Quantum phase transitions (QPTs), or singularities in the evolution of the 
ground state properties of a system as a Hamiltonian parameter is varied, 
have been extensively studied for various many-body systems {e.g., Refs. [1- 
3]). Recently, analogous singular behavior has been noted for states throughout 
the excitation spectrum of certain many-body models [4-10], namely the Lip- 
kin model [11] and the interacting boson model (IBM) for nuclei [12]. These 
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singularities have been loosely described as "excited state quantTim phase 
transitions" (ESQPTs) [9]. In this article, we more closely and systematically 
examine the characteristics of such excited state singularities as phase transi- 
tions, to provide a foundation for future investigations. It is found that excited 
state quantum phase transitions occur in a much broader class of many-body 
models than previously identified. 

Ground state QPTs are characterized by a few distinct but related proper- 
ties. The QPT occurs as a "control parameter" ^, controUing an interaction 
strength in the system's Hamiltonian H{^), is varied, at some critical value 
C = Cc- For specificity, we take the Hamiltonian to have the conventional form 
H{^) = (1 — C,)Hi + ^H2. At the critical value: (1) The ground state energy 
Eq is nonanalytic as a function of the control parameter at ^ = ^c- (2) The 
ground state wave function properties, expressed via "order parameters" such 
as the ground state expectation values {Hi)o or {H2)o, are nonanalytic at 
^ = ^c- These two properties are not independent, since the evolution of the 
ground state energy and that of the order parameters are directly related by 
the Feynman-Hellmann theorem [13], which gives dEo/d^ = {H2)o~ {Hi)o- (3) 
The gap A between the ground state and the first excited state vanishes at 
^ = ^c- (Here we consider only continuous phase transitions. More specifically, 
the systems considered in this article undergo second-order phase transitions, 
in which discontinuity occurs in the second derivative of the ground state en- 
ergy and the first derivatives of the order parameters.) Singularities strictly 
only occur for an infinite number of particles in the many-body system, but 
precursors can be observed even for very modest numbers of particles. For 
finite particle number N, the defining characteristic of the QPT is therefore 
not the presence of a true singularity but rather well-defined scaling behavior 
of the relevant quantities towards their singular large- limits [14]. 

For the systems which exhibit excited state QPTs, the vanishing gap between 
the ground state and first excited state at the ground state QPT does not occur 
in isolation. Rather, there is a bunching of levels near the ground state, that 
is, a vanishing of the average level spacing A or an infinite local level density 
p = A^^. The infinite level density, moreover, propagates to higher excitation 
energy (as illustrated for a two-level fermionic pairing model in Fig. 1) as the 
order parameter is varied from hence the concept of a continuation of the 
QPT to excited states. The singular level density occurs simultaneously with 
singularities in other properties of the excited states \k), such as the level 
energy Ek and the expectation values {Hi)k and {H2)k- 

First, we review the essential properties of the two-level pairing many-body 
models, for both bosonic and fermionic constituents (Sec. 2). We find that 
ESQPTs are universal to these models, suggesting that the ESQPT phenom- 
ena may be broadly relevant, at least to systems dominated by pairing in- 
teractions. The semiclassical analysis of a "sombrero" potential provides a 
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Fig. 1. Excitation energies for the two-level fermionic pairing model (2.7) with par- 
ticle number A'^ = 100, at half filling and zero seniority, as a function of the control 
parameter ^. 

basis for understanding many of the properties of the quantum many-body 
ESQPT [6,9]. The semiclassical analysis of Refs. [6,9] is extended in Sec. 3 
to address several properties relevant to the definition of phase transitions. 
In particular, the singularity in the eigenvalue spectrum and the finite-size 
scaling behavior for the ESQPT are determined at the mean field level. Nu- 
merical calculations for the full quantum problem are considered in Sec. 4, 
where we investigate manifestations of the ESQPT in the excitation spectrum 
and in the properties of "order parameters" for the excited states. Finally, we 
consider the ESQPT as a boundary between quahtatively distinct "phases" 
(Sec. 5). The relationship between the U(n + 1) two-lcvcl boson models and 
the two-level pairing model is established for arbitrary dimension in the ap- 
pendices, where some further mathematical definitions and identities are also 
provided for reference. 



2 Bosonic and fermionic two-level models 

Ground state QPTs have been studied extensively {e.g., Refs. [2,15-17]) for 
the two-level boson models, or s-b models, defined in terms of a singlet bo- 
son s^^^ and a (2L -|- l)-fold degenerate boson 6^^^ [Fig. 2(a)]. Models in this 
class include the U(6) interacting boson model (IBM) for nuclei (L = 2) [12], 
which is defined in terms of s^'^^ and S"^^ bosons, and the U(4) vibron model 
for molecules {L = 1) [18]. Also, the Lipkin model [11] has several isomorphic 
realizations, defined variously in terms of systems of interacting fermions, in- 
teracting spins, or interacting bosons (Schwinger realization). This last real- 
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Fig. 2. Single-particle level degeneracies for the various classes of two-level models 
considered: (a) the s-h boson models, (b) the more general two-level bosonic pairing 
models, and (c) the two-level fermionic pairing models. 

ization falls into the two-level boson model categorization, as the L = case. 
So far, excited state QPTs have been considered in the Lipkin model [5, 6] and 
the IBM [7-9], both of which are examples of s-h two- level models. 

The s-h two-level models are described by the U(n -|- 1) algebraic structure 



U(n + 1) D < 



SO(n+ 1) 
U(n) 



> D SO(n) D S0(3), (2.1) 



where n = 2L+1. The U(n+ 1) generators are given in tensor form by (s^ x sy^^ , 
(s^ X bY^\ {b^ X sY^\ and {b'^ x (see Appendix A for detailed definitions). 
If the Hamiltonian is simply taken as the Casimir operator (A. 11) of either 
of the subalgebras, SO(n -|- 1) or U(n), a dynamical symmetry is obtained. 
The U(n) symmetry is geometrically related to the n-dimensional harmonic 
oscillator, the SO(n-|- 1) symmetry to the n-dimensional rotator- vibrator {e.g., 
Ref. [19]). 

The ground state QPT in the two-level boson models arises as the Hamiltonian 
is varied linearly between the two dynamical symmetries, for instance, by 
varying ^ in the Hamiltonian 

H = ^^AT, - ^{s^b + b^~s) ■ {s^b + b^~s), (2.2) 

where = (-)H^ ■ b is the 6-boson occupancy, TJ^^ = (_)A-m7-W^ g^^^j 
[/(A) . y{A) ^ (-)^(2L + 1)1/2(A X 5)(°). This Hamiltonian yields the U(n) 
symmetry for ^ = and the SO(n + 1) symmetry for = 1. The Hamiltonian is 
invariant under the common SO(n) algebra in (2.1) and therefore conserves a 
(2L-|- l)-dimensional angular momentum quantum number v. As ^ is increased 
from ^ = 0, the increasing strength of the interaction between s and b levels 
changes the structure of the ground state from a pure s-boson condensate to a 
condensate involving both types of bosons. For asymptotically large values of 
the total particle number N = Ns + N/,, the change is abrupt. A second-order 
ground state QPT is well known to occur for ^c= 1/5, with all the properties 
enumerated in Sec. 1. [The coefficients in (2.2) are scaled by appropriate pow- 
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ers of to guarantee that the location of the critical point is independent of 
in the large limit.] With more complex interactions in the Hamiltonian, 
first-order QPTs, such as the physically important U(5)-SU(3) phase transi- 
tion in the IBM, may also be obtained [2, 15, 20]. The conditions under which 
such first-order phase transitions occur in an arbitrary \J{n + 1) model are 
outlined in Ref. [21]. However, only second-order ground state QPTs will be 
considered here. 

We observe, moreover, that the U(n + 1) two-level boson models are special 
cases of an even larger family of models, the two-level pairing models with 
quasispin Hamiltonians. Two-level pairing models can be defined for systems 
of either bosons or fermions. The two-level pairing models undergo a second- 
order ground state QPT [22]. Therefore, it is natural to consider the possibility 
that excited state QPTs may occur within the context of this broader family 
of models as well. 

The quasispin pairing Hamiltonian is of the form 

(E c]mCjm) + ^ E (E 4^' 4m' ) (E CimCjm) , (2.3) 

j ^ m j'j m' ^ ^ m ^ 

where the summation indices j and j' run over the single-particle levels, and 
m and vn! run over their substates. The Cj may represent either bosonic op- 
erators ft^^^'* and h^2^^ [Fig- 2(b)] or fermionic operators operators ci^l^^ and 
02^^ [Fig. 2(c)], as appropriate. Although the Hamiltonian (2.3) superficially 
appears quite different from the U(n -|- 1) two-level boson model Hamilto- 
nian (2.2), the two are in fact equivalent [23,24]. The detailed relationship 
between the models is established for arbitrary n in Appendix A. 

It is well known that the pairing Hamiltonian (2.3) can be expressed in terms 
of the generators and Sjz of a quasispin algebra (A. 2), as 

E = E^.(2'5,. T %) ^Y.Gy,Sy^S,., (2.4) 
J j'j 

where f2j = (2j + l)/2, and the upper and lower signs apply in the bosonic and 
fermionic cases, respectively. The algebra is either an SU(1, 1) algebra if the op- 
erators are bosonic [25] or an SU(2) algebra if the operators are fermionic [26]. 
However, the pairing models are also characterized by an overlaid U(ni -|- n-i) 
algebraic structure, described further in Ref. [27], either 

D SOi(ni) ® S02(n2) D SOi2(3) (2.5) 



U(ni -t- n2) D < 



SO(ni + na) 
Ui(ni) ® U2(n2) 
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in the bosonic case (with rii — 2Li + 1 and n2 — 2L2 + 1) or 

D Spi(ni) Sp2(n2) D SUi2(2) (2.6) 

in the fermionic case (with rii = 2ji + 1 and 77,2 = 2j2 + 1), directly generahzing 
the U(n + 1) algebraic structure (2.1) of the s-b boson models. The U(ni + 77,2) 
generators arc of the form (c| x CiY'^\ (c| x C2Y'^K {cl x CiY^\ and (c^ x c^Y^^- 
The SOi(ni) and 802(^2) [or Sp]^(ni) and Sp2(n2)] algebras provide conserved 
ni-dimensional and n2-dimensional angular momentum quantum numbers {vi 
and V2), which are equal to the seniority quantum numbers defined in the 
quasispin formulation. 

The ground state QPT in the general two-level pairing models is between the 
SO(ni + 112) or Sp(ni + 712) dynamical symmetry and the Ui(rii) ® 1)2(^2) 
dynamical symmetry. To choose a transitional Hamiltonian for the general 
pairing models consistent with the Hamiltonian already used for the s-b bo- 
son models, we observe that the Hamiltonian (2.2) may be reexpressed (see 
Appendix A) in pairing form as 

H = ^^^2 + ^(-)^+^(^i+ ± S2^m- ± ^2-), (2.7) 

to within an additive constant, where the full relation is given explicitly 
in (A. 17). With this form of Hamiltonian for the pairing models, the ground 
state QPT again occurs at = 1/5. 

Since QPTs occur in the limit of large particle number, an important distinc- 
tion arises between bosonic and fermionic models. Arbitrarily large particle 
number can be achieved in the bosonic models, even for fixed level degenera- 
cies, simply by increasing the total occupancy. For a fermionic model, how- 
ever, the total occupancy is hmited by Pauh exclusion to the total degeneracy 
[(2ji + 1) + (2j2 + 1)] ■ Therefore, the hmit of large particle number can only be 
achieved if the number of available substates in each level is simultaneously 
increased. For two fermionic levels of equal degeneracy (ji —j2 =j), half- filling 
is achieved for iV = 2j -|- 1. 

It will be convenient to make extensive use of the U(3) two-dimensional vibron 
model [28-30] for illustration in this article. The U(3) vibron model is the 
simplest two-level model which still retains a nontrivial angular momentum 
or seniority quantum number (unlike the Lipkin model). 

First, note that single-particle levels in bosonic pairing models [Fig. 2(a,b)] are 
only restricted to odd degeneracies {i.e., 2L + 1 with L integer) if a physical 
three-dimensional angular momentum subalgebra S0(3) is required in (2.1) 
or (2.5). The pairing interaction only requires the definition of time-reversed 



\J{ni + 772) D 



Sp(72i + 712) 
Ul(7li) (g) U2(7l2) 
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Fig. 3. Eigenvalue spectra for the U(3) vibron model I = states (A'^ = 100), for 
several specific values of the Hamiltonian parameter ^. Eigenvalues are plotted with 
respect to the scaled excitation quantum number k/N. 

pairs. It therefore suffices to have an "M" quantum number, with pairs ±M, 
without necessity for an "L" quantum number. The pairing interaction can 
therefore be defined for an even number of bosons, and the interaction within 
each level is described by SO(n) with n even. ^ Bosonic levels of even degen- 
eracy arise naturally in problems lacking three-dimensional rotational invari- 
ance. 

The U(3) vibron model may be obtained by considering the U(4) vibron model 
(L — l) and eliminating the substate 6o^\ This leaves a U(3) algebraic struc- 
ture, with S0(3) and U(2) dynamical symmetries. The geometrical coordinates 
associated with the U(4) model describe three-dimensional dipole motion (as 
in a linear dipole molecule). However, eliminaton of ^o^"* "freezes out" motion 
in the z direction, so the U(3) model instead describes two-dimensional motion 
in the xy plane. The U(2)-SO(3) transitional Hamiltonian, in Casimir form, 
is [28,29], 

'MD+D. + D.D+)+f], (2.8) 



N ^ m 



where D± = ±\/2(6!lj_iSo — So^ipi)- This is the two-dimensional equivalent 



^ However, pairing is not well-defined for the converse situation, a fermionic level of 
odd degeneracy. With an odd number of substates, one ("m = 0") must necessarily 
be its own conjugate under time reversal. Creation of a time-reversed pair involving 
this substate is Pauli forbidden. The corresponding algebra, Sp(n) with n odd, is 
not defined. 
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of (2.2), to within an additive constant. The conserved two-dimensional angu- 
lar momentum is / = fc^i^+i — The eigenvalue spectra of / = states, 
for various values of ^, are shown in Fig. 3. The spectra for the U(2) dy- 
namical symmetry = 0) and the S0(3) dynamical symmetry = 1) have 
simple analytic forms [28] . Note also the spectrum for the ground state QPT 
(e = 0.2). 



3 Semicleissical dynamics 

3.1 Coordinate Hamiltonian 

Each of the many-body models considered in Sec. 2 has an associated classi- 
cal Hamiltonian, defined with respect to classical coordinates and momenta, 
which is obtained through the use of coherent states [2,31,32]. The basic 
properties of the excited state quantum phase transition follow from the semi- 
classical analysis of a double-well potential with a parabolic barrier [Fig. 4(c)] 
or, in higher dimensions, a sombrero potential (also known as the "cham- 
pagne bottle" potential [33]). The semiclassical dynamics for these potentials 
has been studied in depth [33-37], and the connection with ESQPT phenom- 
ena in the Lipkin model and higher-dimensional s-b boson models has been 
made in Refs. [4, 6, 8, 9]. In particular, at the energy of the top of the barrier, 
the classical action undergoes a logarithmic singularity, which leads semiclas- 
sically to the prediction of an infinite level density. Here we do not attempt a 
comprehensive recapitulation of the existing analysis but rather briefly sum- 
marize the essential points and derive some results speciflcally relevant to the 
observables of interest in phase transitional phenomena. 

For the quasispin models of Sec. 2, the two superposed algebraic structures 
(quasispin and unitary) give rise to two alternative sets of coherent states and 
therefore to two realizations of the classical dynamics. The SU(1, 1) or SU(2) 
quasispin algebra yields a one-dimensional dynamics (the phase space is a 
Bloch sphere or hyperboloid [31, Ch. 6]) which is common to all the quasispin 
models. The dynamics arising from the quasispin algebra therefore highlights 
aspects universal to these models, yielding the basic double-well potential 
[Fig. 4(c)] and therefore indicating that all should exhibit an ESQPT at the 
energy of the top of the barrier. In contrast, the coherent states obtained from 
the unitary U(ni -|- ^2) algebra yield a much richer classical dynamics, in nin2 
dimensions, associated with the coset space U(ni -|- n2)/[U(ni) (8) U(n2)] [31, 
Ch. 9] . This more complete dynamics, so far only fully investigated for the s-b 
models [38,39], yields a much more detailed description of the system. The 
dynamics obtained from the quasispin algebra is essentially a one-dimensional 
projection or "shadow" of the full dynamics arising from the unitary algebra. 
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as described by Feng, Gilmore, and Deans [2] for the IBM. In particular, 
the presence of angular degrees of freedom and conserved angular momentum 
quantum numbers have significant consequences for the ESQPT [8,9]. 

First, let us summarize the classical Hamiltonian obtained from the \J{n + 
1) coherent states for the s-b model. The classical Hamiltonian acts on n 
coordinates and their conjugate momenta. However, for the SO(n)-invariant 
interaction in (2.2), the Hamiltonian is invariant under rotations in the n- 
dimensional space and can therefore be expressed solely in terms of a radial 
coordinate r, its conjugate momentum pr, and a conserved angular kinetic 
energy T^{v), as [9, 17,39] 

H-^-^\pI + r-'Uv)] + ^^[tY + Uv)] + ^-^r' + ir\ (3.1) 

where T^{v) has eigenvalue f (f + n — 2) and the coordinate r is defined only 
on the domain < r < \/2. ^ The eigenvalue problem for (3.1) therefore has 
the form of a radial Schrodinger equation with a quadratic-quartic potential, 
except for the appearance of the position-dependent kinetic energy term pro- 
portional to r^p^. For the one-dimensional case, i.e., the Lipkin model, the 
centrifugal term is not present, and the coordinate and momentum are more 
aptly denoted by x and p, so 

" ~2N^^ + + ^ ' (3.2) 

where here both negative and positive values of the coordinate x are allowed 
(-V2<x<+\/2). 

The role of h'^/{2m) in the usual Schrodinger equation is taken on by the 
coefficient of or in (3.1) or (3.2). We therefore make the identification 
N~^, with the coordinate-dependent mass m(a;) = (1 — ^ -|- 2^a;^)~^. 

The forms assumed by the quadratic-quartic potential in (3.1) or (3.2), V{x) ~ 
(1 — 5^)x'^/2 + ^x^, are summarized for convenience in Fig. 4(a-c). For the 
radial problem, of course, only the positive abscissa is relevant. For ^ < 1/5, 
the potential has a single minimum, at a; = 0, which is locally quadratic. For 
^ = 1/5, the critical value for the ground state QPT, the potential is pure 
quartic. For ^ > 1/5, the familiar double-well potential is obtained (or the 
sombrero potential for n> 1). For the Hamiltonian (3.1) or (3.2), the zero in 
energy is such that the top of the barrier is always at £^ = 0, independent of ^. 



^ In obtaining (3.1) from Ref. [39], a scaling transformation r — > N^/^r has been 
made, and a constant term of order 1/N has been suppressed. 
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3.2 Singular properties of the action 

The main semiclassical features of levels at energies near the top of the bar- 
rier are obtained by noting that for E = the classical velocity v{E,x) — 
[2[E — V{x)]/m{x)]^^'^ locally vanishes at the top of the barrier (x = 0). While 
indeed the classical velocity also vanishes at the ordinary linear turning points 
of a potential well, the vanishing slope at the top of the barrier presents a qual- 
itatively broader "flat" region over which the classical velocity is small. Thus, 
the semiclassical motion has a long "dwell time" in the vicinity of a; = 0. This 
leads to two essential results, namely (1) an infinite period r = § v{E,x)^^ dx 
for classical motion across the top of the barrier and (2) strong localization 
of the semiclassical probability density P{x) oc v^E^^x)^^ at the top of the 
barrier [6]. 

The first-order semiclassical analysis provides a simple guidemap to the prop- 
erties of the spectrum as a whole and also provides an explanation for the 
singularity in level density as the top of the barrier is approached. We con- 
sider the one-dimensional problem (3.2), but the results apply equally to the 
radial problem (3.1) with f = 0. For the Hamiltonian (3.2), the usual first-order 
WKB quantization condition [40] becomes 



with A; = 0, 1, . . ., where the action S = §pdx over a full classical period of 
motion is given by the integral 



between classical turning points xi{E) and X2{E).^ The action depends upon 
^ variously through m(x), l^(x), and the turning points. 

The quantization condition (3.3) implicitly gives the adiabatic evolution of 
the energy Ek{C) of ^ given level with respect to the parameter ^. Since (3.3) 
enforces that S'(^; E) be constant if k is held fixed, the curve describing Ek{^) 
is simply a contour of S{S,\ E) in the ^-E plane. These contours, calculated 
numerically for the Hamiltonian (3.2) [or (3.1) with v = 0] are plotted in Fig. 4. 

^ Some bookkeeping issues naturally must be taken into account in the one- 
dimensional double- well problem [Fig. 4(c)]. For E < 0, i.e., below the barrier, 
the two wells are classically isolated. Applying the quantization condition with S 
evaluated over one of the wells in isolation is equivalent to counting only states 
of one parity (symmetric or antisymmetric). For E > 0, applying the quantization 
condition with S evaluated over the full well counts states of both parity. Questions 
as to the proper transition between the regimes E < and E > are somewhat 
artificial, since the validity conditions for (3.3) break down at E'RiO. 



SiC;E)^ik + l)2nN 



-1 



(3.3) 




(3.4) 
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Fig. 4. Contour plot showing the global structure of the classical action S{(,; E) for 
the geometric Hamiltonian (3.1) or (3.2), through the different regimes determined 
by the shape of the quadratic-quartic potential energy function (2), which is shown 
for (a) ^ < Co (b) ^ = Cc) and (c) ^ > Cc- The individual contours are related 
semiclassically to the evolution of the level eigenvalues Ek{^). 

A compression of energy levels at = is visible qualitatively even here. [In 
Fig. 1, the Ek{C,) are plotted as excitation energies and therefore cannot be 
compared directly with Fig. 4. More appropriate plots for comparison may be 
found in the following section, e.g., Fig. 7(a).] The derivative dEk/d^ along a 
single contour of S{^;E) is plotted in Fig. 5(a). Note that dE^/d^ undergoes 
a singularity in which dE^/d^ — > but d'^E^/d^'^ — > ±oo, at a critical value 

In semiclassical analysis, the gap or level density is directly related to the 
classical period. Prom the quantization condition (3.3), it follows that the 
semiclassical estimate of the gap between adjacent levels (A = dE^/dk) is 
A{E) = 27rN-\dS/dE)-\ By differentiation of (3.4), the gap is simply 
A(^) = 27rN-^T-^. As already noted for the ESQPT [6], the period r be- 
comes infinite at E'^O and, equivalently, the gap A(£') vanishes. An exphcit 
calculation of {dS/dE)~^ as a function of E for the classical Hamiltonian (3.2) 
is shown in Fig. 5(b). Note that {dS/dE)~^ undergoes a singularity in which 
(dS/dE)-^ but (d'^S/dE^)-^ ±oo, at the critical energy E^ = 0. 

For nonzero angular momentum v, the origin (r = 0) is classically forbidden 
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Fig. 5. Singularities in derivatives of the classical action (3.4) for the geometric 
Hamiltonian (3.1) or (3.2). (a) The derivative dE/d£, along a countour of S{^;E) 
(Fig. 4) , related semiclassically to the adiabatic evolution of the level energy, (b) The 
inverse of the partial derivative (dS/dE)^, proportional to the semiclassical estimate 
for the gap. 

due to the centrifugal term in (3.1), which causes the wave function probability 

near the origin to be suppressed. This mitigates the effects just described, by 
masking the top of the barrier and precluding the long semiclassical dwell time 
at the origin [9]. The dependence of the Hamiltonian (3.1) on v is through the 
coefficient of the centrifugal term, which is proportional to T^{v) /N'^ ~ (v/N)^. 
Therefore, the phenomena associated with the ESQPT can be expected to 
be suppressed for sufficiently large v at any given value of N. On the other 
hand, the angular momentum effects at any given value of v are negligible 
for sufficiently large N. That is, the signatures of the ESQPT persist for 
small V {v/N « 1) and only disappear for v/N of order unity (as illustrated 
quantitatively in Sec. 4.1). 

3.3 Asymptotic spectrum 

Let us now consider more precisely the form of the singularity in the spectrum 
in the immediate neighborhood of the ESQPT. As the wave function becomes 
increasingly well- localized near the top of the barrier for £■ — > 0, it should 
become an increasingly good approximation to treat the barrier as a pure 
inverted oscillator potential, V{x) = —Ax^. The position-dependent kinetic 
energy term (oc x^p^) also becomes irrelevant. 

In the action integral (3.4), the classical turning point at the barrier is xi{E) ~ 
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[E /Ay/"^ for < 0, or for E >0 integration simply extends to the origin. The 
distant turning point X2{E) is a slowly varying function of E which does not 
contribute to the singularity, so we may take it to be a constant. (In any case, 
for the actual potential, the approximation of a pure parabolic barrier breaks 
down well before the distant turning point is reached.) The action integral for 
£■ > is therefore 



2 

h Jo 

4 i-{A/E)y^X2 



h-^S{E) = =: r dx [2m{E + Ax^)Y'^ 
n Jo 



—E [ 
fko Jo 



(3.5) 



where, for the inverted oscillator Hamiltonian H — [H^ / (2m)]p^ — Ax^ , we have 
defined hu; — 2[/i^/(2m)]^/^74^/^ by analogy with the conventional harmonic 
oscillator. 

Expanding this action integral [41, (2.271.3)] for E'faO yields 

h-^S{E) = (-E \ogE + ao + aE + ■■■), (3.6) 

where aQ and a are constants, i.e., depend only on the potential parameters 
A and X2. An essentially identical result is obtained for E < 0, with the re- 
placement > l^^l [41, (1.646.2)]. The singular behavior for energies near the 
top of the barrier therefore arises from the E log E term. ^ The quantization 
condition (3.3) takes on the form 

-ElogE + aE ^ ^27rhu}{k- kc), (3.7) 

where = is obtained for k = k^. If the energy dependence in (3.7) is trun- 
cated at the terms shown, i.e., linear order in E, this quantization condition 
can be solved for E{k) in terms of the Lambert W function, by (B.5), yielding 
E{k) — —2'Khw{k — kc) /W[—e~°'2'Khuu{k — k^)] - The relevant properties of the 
W function are summarized in Appendix B. 

For the Hamiltonian (3.2), the top of the barrier is described by an oscillator 
constant which may be read off from the coefficients of the and x^, giving 

Sf£)V2 

hu; = (3.8) 

where 

5(0^(1-0(1-50- (3.9) 



Since the Schrodinger equation for a pure parabolic barrier is exactly solvable 
in terms of parabolic cylinder functions [42], the \E\ log|£^| dependence can also be 
obtained by explicitly matching this solution for the wave function in the vicinity 
of the barrier to asymptotic WKB wave functions away from the barrier [43] . 
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The oscillator constant thus depends upon both ^ and A^. The same function 
interestingly, also enters into the ground state QPT scaling properties, as 
obtained by the continuous unitary transform method in Ref. [44]. The semi- 
classical estimate for the eigenvalue spectrum in the vicinity of the ESQPT is 
therefore 

where a will contain a dependence on ^ as well. Differentiation with respect 
to k, making use of (B.4), yields a semiclassical estimate 

^^"^'^'^^ = -W[-e-^2nEiO^/^ik-k.)/N] + l ^^'''^ 
for the energy gap between adjacent excited states. 

Since the excitation quantum number and particle number enter into the quan- 
tization condition (3.3) together in the combination k/N, the spectrum and 
finite-size scaling properties are inextricably linked at the semiclassical level. ^ 
The expression (3.11), considered as a function of at fixed k, provides an 
estimate for the scaling of the gap at the {k — A;c)-th eigenvalue above or below 
E = 0. The large-iV behavior follows from the known asymptotic form (B.2) 
of the W function as a sum of logarithms for x — > 0~ (see Fig. B.l). The 
values of x relevant to (3.11) in the vicinity of the ESQPT are of the order 
X ~ — A^~^. The asymptotic form (B.2) provides a good approximation to 
W{x) for reasonable N, e.g., accurate to 1% by 10^. 

For very large A^, the scaling behavior is in principle even simpler. The \og{—x) 
term in (B.2) outgrows the log[— log(— x)] term as a; — > 0~. With this logarith- 
mic approximation, an extreme asymptotic estimate 

7VA~ '-^^ -H^Htt)!^! (312) 

log(A; - k^) - log + log[27rS(0 V2] -a + 1 log A^ ^ ' ^ 

is obtained, recovering the logarithmic scaling noted by Leyvraz and Heiss [6]. 
However, even for N ~ 10^°, the approximation W{x) ~ log(— x) yields an 
error of > 10% and therefore is of limited quantitative value for systems of 
typical "mesoscopic" size. 

Note that the quantization condition as given in (3.3) is derived under the 
assumption that the classical turning points are well separated (by several 
de Broglie wavelengths) and that the potential is locally linear at these turn- 
ing points [40]. This suffices for the analysis of levels which are not close in 



^ For the ground state QPT, the semiclassical potential is quartic [Fig. 4(b)]. A 
simple application of the WKB formula gives a dependence E{k/N) ~ (k/N)^/^, 
which simultanously yields both the spectrum Ek ~ k*^^ [Fig. 3 (.^ = 0.2)] and the 
scaling A ~A^-4/3 (Sec. 4.2). 
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energy to the top of the barrier. However, for E 0, the barrier presents a 
quadrafic classical turning point. (Equivalently, the linear turning points on 
either side of the barrier approach each other, violating the assumption of 
sufficient separation.) For accurate quantitative analysis of the levels imme- 
diately surrounding E — Q, the more general phase- integral method must be 
applied [45]. For a smooth, symmetric doTiblc-wcll potential, the phase-integral 
method yields an approximate quantization condition [45, (3.47.1)] 

h~^S{E) = 27r{k + I) - + 2/?^' ± arctanexp(-ir), (3.13) 

with k an integer, where the various phases appearing on the right hand side 
are defined in Ref. [45]. The full derivation involves the evaluation of contour 
integrals on the complex extension of the coordinate axis and the consideration 
of complex- valued turning points for energies just above the barrier [45] . Quan- 
titative solution of the problem is considered in detail in Refs. [35,37,46,47]. 
The effects of these corrections (3.13) relative to (3.3) are explored in Ref. [33]. 
The corrections are essential to the treatment of the first few eigenvalues above 
or below the barrier. However, here we are instead interested in extracting the 
basic nature of the singularity from the dependence of S{E) on E in the 
vicinity of £' = 0, for which the simple quantization condition (3.3) suffices. 



4 Quantum properties 

4-1 Eigenvalue spectrum 

In a ground state QPT, the singular behavior of the system is simultaneously 
reflected in the eigenvalue spectrum (ground state energy and gap) and in 
the order parameters. From the preceding semiclassical analysis (Sec. 3), it is 
to be expected that a similar variety of interconnected phenomena occur at 
the ESQPT, and this is indeed borne out by the quantum calculations. Of 
course, the analogy between ground state QPT and ESQPT is far from exact, 
so let us now examine the results for spectra and order parameters obtained 
numerically from the full quantum calculation, to elucidate both the analogy 
with the ground state QPT and the apphcability of the semiclassical results 
of Sec. 3. 

While the ground state QPT may only be traversed by varying a Hamiltonian 
parameter, the locus of the ESQPT is a curve in the two-parameter space 
defined by the Hamiltonian parameter and the excitation energy (as along 
the dense band in Fig. 1). Therefore, the ESQPT may be crossed cither "hor- 
izontally", by varying ^, or "vertically", by varying the excitation quantum 
number k (or, equivalently, the energy E) of the level being examined. 
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The energy spectrum consists of the set of eigenvalues £'7v,^,fc,A, which contain 
dependences on several quantities: the system size A^, the Hamiltonian pa- 
rameter ^, the excitation quantum number k, and other conserved quantum 
numbers A (angular momenta or seniorities in the present models). For large 
N, however, k/N and A/N become essentially continuous variables. In the 
preceding section, it was seen that semiclassically the energy depends upon 
the quantum numbers only through these combinations k/N and A/N. We 
are therefore largely interested in the properties of the spectrum given by the 
function E{^, k/N, A/N) of three quasi-continuous variables [17]. The depen- 
dence of the spectrum on interaction, excitation quantum number, and angu- 
lar momentum is contained in the dependence of E[^, k/N, A/N) on its three 
arguments.^ Furthermore, note that the dependence on the argument k/N 
implicitly contains information not only on the excitation spectrum [when the 
function is considered as E{k) at fixed N] but also on the finite-size scaling be- 
havior [when the function is considered as E{l/N) at fixed k]. The properties 
of E{^, k/N, A/N) in the vicinity of the ground state, that is, for k/N <^ 1, 
have been studied in detail, at least for the s-b models. Here, instead, we are 
considering the regime k/N r^l. 

First, let us establish the common ground between the various models under 
consideration (Sec. 2), by a simple comparison of the energy spectra. Calcu- 
lations are shown in Fig. 6 for the Lipkin model [Fig. 6(a)], the U(3) vibron 
model [Fig. 6(b)], a bosonic pairing model with equal degeneracies for both 
levels {Li = L2 = I) [Fig. 6(c)], and a fermionic pairing model with equal 
degeneracies (ji = j2 — 9/2) [Fig. 6(d)]. The calculations are all for a fixed, 
modest particle number (N — 10), so that individual eigenvalues are clearly 
distinguishable. In the comparison, we must distinguish the invariant sub- 
spaces of states for each model. Each eigenstate of the Lipkin model contains 
only even-Nh or odd-A^f, components and is thus characterized by a grading 
quantum number g with values and 1 (g^ Ni, mod 2) or, equivalently, the 
parity tt = (— )^. The vibron model states are characterized by the angular mo- 
mentum / = 0, ±1, . . . , ±A^. The bosonic and fermionic pairing model states 
are characterized by seniority quantum numbers for each single-particle level, 
namely, vi and V2- 

Note the essentially identical evolution, with respect to ^, of the even-parity 
{g — 0) states of the Lipkin model, the zero angular momentum {I — 0) states 



^ In the Hamiltonians (2.2) and (2.7), the coefficients of the one-body operators 
are scaled by N and the coefficients of the two-body operators are scaled by N'^. 
Often a Hamiltonian normalization differing by an overall factor of N is instead 
used, e.g., H = {I - ^)Nh - {C/N){s^b + bU) ■ (stfo + b^s) for the s-b model. For 
the normalization (2.2) or (2.7), -Eiv.^.fc.A docs indeed approach a limiting value 
E{^, k/N, A/N) as N^oo, by (3.3). However, for the alternate normalization it is 
actually iV~^£^7v,^,fc,A which approaches a limiting value as N^oo. 
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Fig. 6. Eigenvalues for (a) the Lipkin model (Schwinger realization), (b) the U(3) vi- 
bron model, (c) the bosonic pairing model (Li = L2 = 1), and (d) the fermionic pair- 
ing model (ji = j2 =9/2), as functions of the coupling parameter ^, all for total par- 
ticle number N = 10. For the Lipkin model, both even-parity (solid curves) and odd- 
parity (dashed curves) levels are shown. For the other models, only the lowest angu- 
lar momenta or seniorities are shown. A diagonal contribution ^(N -\-2Li + 2L2)/N 
has been subtracted from the Hamiltonian (2.7) for the bosonic pairing model [27]. 

of the vibron model, and the zero seniority [{V1V2) — (00)] states of both the 
bosonic and fermionic pairing models (solid curves in Fig. 6). The ground state 
energy is near constant, with EQ^{)^ioi ^<^c and decreases to < — 1 for ,^ = 1. 
The highest eigenvalue decreases approximately linearly with ^, from 1 to 0. 
Various quahtative features associated with the ESQPT occur at £■ for 
^ > for these models. Note especially the inflection points for these levels 
(solid curves) as well as the change in the pattern of degeneracies between 
different seniorities (or parities or angular momenta) at E^iQ. 

The major differences among the models lie in the degeneracy patterns at 
nonzero seniority, which depend upon the speciflc algebraic properties of the 

individual models [27]. At present, we will limit consideration of angular mo- 
mentum effects to the s-h models, since for these only one angular momentum 
quantum number is involved, and / in the U(3) vibron model (Sec. 2) serves 
as a natural example for illustration. 
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Fig. 7. Angular momentum dependence of spectral properties for the U(3) vibron 
model {N = 100). (a,b) Evolution of eigenvalues with ^ for / = and I = 25, i.e., 
l/N = 1/4. (c) Dependence of the gap on excitation energy, as in Fig. 8(b), for 
various / (0</<25). 

The semiclassical analysis of Sec. 3 provided a simple set of predictions (Fig. 5) 
for the singular behavior of E{^,k/N,A/N) as the ESQPT is crossed both 
"horizontally" ^-iid "vertically" [E{k/N)]. Namely, E{^) undergoes a 

singularity in which the slope sharply approaches zero {dE/d^ 0) [Fig. 5(a)] 
but with a curvature which becomes infinite and reverses sign {d'^E/d^^ — >• 
±oo), yielding a special divergent form of inflection point, as ^ — > A similar 
singularity is expected in E{k/N) [Fig. 5(b)] at the critical energy. 

The actual diagonalization results at finite N show clear precursors of this 
form of singularity in £■ as ^ is varied. Even for the small system size {N — 10) 
considered in Fig. 6, each eigenvalue E{^) undergoes an inflection [Fig. 6 
(solid curves)] at an energy close to the expected critical energy, i.e., E^ — Q 
for the Hamiltonians used. The derivative dE/dC, is shown for larger boson 
number [N = 100 and 1000) in Fig. 8(a), for the U(3) vibron model / = 
states. The second derivative d'^E/d^'^ is also shown (inset). The expected dip 
dE/d^ — > and divergent inflection d^E/d^"^ — > ±oo both are present and 
become gradually sharper with increasing N. 

For nonzero / in Fig. 6(b), the inflection points in the eigenvalues as functions 
of ^ are washed out, as expected from the semiclassical analysis (for l>0 the 
centrifugal term suppresses the probability density near r — 0, mitigating the 
effect of the barrier). For A^ = 10, the inflection points disappear even for the 
very lowest nonzero / values [dashed curves in Fig. 6(b)]. Also for = 10, 
the inflection points are suppressed for the negative parity {g = 1) states of 
the Lipkin model in Fig. 6(a). Here a similar mechanism applies: negative 
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Fig. 8. Evolution of excited level energies and the order parameter (iVf,) across the 
ESQPT, as traversed both by varying ^ (left) and by varying E (right), i.e., "hori- 
zontally" and "vertically" . Calculations are shown for the U(3) vibron model / = 
states, with N = 100 (dashed curves) and 1000 (solid curves), (a) The derivatives 
dE/d^ and d'^Ejd^^ (inset), for a specific excited level (A;/A^ = 0.2). (b) The deriva- 
tive dE/d{k/N) or, equivalently, the scaled gap A^A, and d'^E/d{k/N)'^ (inset), for 
^ = 0.5. (c) The order parameter {N},) (rescaled by N) as a function of ^ for the same 
level as in panel (a), (d) The order parameter (A^b) (rescaled by A'') as a function of 
excitation energy, for the same ^ value as in panel (b). The discrete eigenstates are 
resolved at the expanded scale shown in the inset. 

parity states posess a node at a; = 0, and the effect of the parabolic barrier 
at a; = is therefore again reduced. (To this extent, the grade in the Lipkin 
model is a surrogate for the angular momentum in the higher-dimensional 
boson models. The formal relation is given in Appendix A.) Compare also the 
curves for nonzero seniorities in Fig. 6(c,d). While the change in behavior of 
the eigenvalues between / = and nonzero / seems to be rather abrupt for the 
= 10 illustration, it must be borne in mind that the relevant parameter for 
the semiclassical description was noted to be l/N, which can only be varied 
very coarsely when A^ = 10. The more gradual evolution of the ESQPT with 
l/N, as obtained for larger A^, is considered further below. 

The properties of the spectrum as the ESQPT is traversed "vertically" by 
varying the excitation quantum number for a single fixed Hamiltonian pa- 
rameter value ^ are explored in Fig. 8(b), again for the U(3) vibron model 
with A^ = 100 and 1000, now at the specific parameter value ^ = 0.5. The 
singularity in E{C,,k/N,A/N) with respect to k/N gives rise to the origi- 
nal, defining property of the ESQPT, namely the vanishing gap or infinite 
level density. The gap is simply the change in energy for a unit change in k 
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quantum number, so in the limit where k/N is taken as a quasi-continuous 
variable we have dE{^,k/N, A/N)/d{k/N) = NA{^, k/N, A/N). The gap is 
shown as a function of energy, rather than of k, in Fig. 8(b), so that the en- 
ergy in the spectrum at which the precursors of the singularity occurs can 
be compared with the expected critical energy Ec — 0. The second derivative 
d'^E/d{k/Ny is also shown (inset). The qualitative features dE/d{k/N) — > 
and d'^E/d{k/N)^ ±00 expected from the semiclassical analysis are indeed 
realized, more sharply with increasing N. 

The inflection point of E with respect to A; at £■ = (though not its singular 
nature) is also immediately visible simply by inspection of the / = spectra 
obtained for various ^ (Fig. 3). The spectra are concave downward with respect 
to k below E = and concave upward above this energy. At the SO (3) limit, 
the entire spectrum falls below E — and constant downward concavity follows 
from the exact formula [28] for the eigenvalues, quadratic in k. Although here 
wc arc considering the dip in dE/dk as a property of the ESQPT in a many- 
body interacting boson model, it should be noted that the dip arising for the 
associated two-dimensional Schrodinger equation is well known as the "Dixon 
dip" [48], with applications to molecular spectroscopy (see also Ref. [30]). 

For nonzero /, as noted above, the relevant parameter governing the disap- 
pearance of the ESQPT is expected to be l/N. The eigenvalue spectrum for 
the U(3) vibron model with N = 100 indeed shows compression of the level 
density at the critical energy for Z = [Fig. 7(a)] and, conversely, no apparent 
compression of level density for large l/N [Fig. 7(b)], where Z = 25 or l/N — l/A 
is shown in this example. (See Ref. [7] for analogous plots for the IBM.) How- 
ever, the gradual nature of the evolution with l/N is seen by considering the 
dip in dE/dk, which becomes continuously less deep and less sharp as l/N is 
increased [Fig. 7(c)]. 



4-2 Finite-size scaling 

The spectroscopic hallmark of the critical point of a QPT is not a vanishing 
gap per se, since the gap never strictly vanishes for finite system size, but 
rather the nature of its approach to zero as increases. It is therefore essen- 
tial to characterize the finite size scaling behavior of the gap in the vicinity of 
the ESQPT. With the Hamiltonian normalization of (2.2), the gap A every- 
where approaches zero with increasing N, so we are actually, more precisely, 
interested in the scahng of the gap at the ESQPT relative to the scahng else- 
where in the spectrum. For states well-separated from both the ground state 
QPT and ESQPT, the scaling is as A ~ A^~^. For states in the vicinity of 
the ground state QPT, the gap vanishes more quickly than A^~^, as the power 
law A ~ A^~^/^. This has been established both numerically and analytically 
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for the various models under consideration [30,44,49-52].^ The gap at the 
excited state QPT also approaches zero more rapidly than A ~ A^~^. This is 
apparent even from the simple plot Fig. 8(b), where A^A is essentially inde- 
pendent of away from the critical energy (compare the curves for = 100 
and N — 1000) but approaches zero with increasing N at the critical energy. 

Let us now examine finite-size scaling more carefully, in particular, to see the 
extent to which the semiclassical expression (3.11) reproduces the scaling be- 
havior. It is not a priori obvious that the semiclassical result (3.11) should 
yield the proper scaling properties for the eigenvalues in the vicinity of the 
ESQPT. Even in the solution of the ordinary Schrodinger equation, the semi- 
classical analysis becomes unreliable for the first few eigenvalues in the vicinity 
of the top of the barrier [33-35] . 

First, in Fig. 9(a), the actual form of the spectrum in the vicinity of the 
ESQPT, obtained by numerical diagonalization, is compared with the semi- 
classical estimate (3.10). Eigenvalues are shown for N = 100 and = 1000. 
Note that kc is simply determined as the value of k for which the energy eigen- 
values cross zero. This must be interpolated between discrete eigenvalues, so kc 
is in general noninteger. The singular logarithmic term in (3.6) has a coefficient 
which is predicted unambiguously from the value of huj (3.8) for the inverted 
oscillator, but no attempt is made here to directly calculate the coefficient a 
of the nonsingular linear term. Rather, a is simply chosen to numerically re- 
produce the linear trend in the eigenvalues in the vicinity of the ESQPT. The 
a value obtained from a limited number of eigenvalues around E — there- 
fore depends somewhat on both A^ and the number of eigenvalues considered. 
The gap, that is the first difference of the eigenvalues in Fig. 9(a), is plotted 
in Fig. 9(b), together with the semiclassical estimate (3.11). The form of the 
singularity is well matched by the semiclassical estimate. (The parameter a 
essentially determines the normalization of the curve A[{k — kc)/N].) The most 
significant deviation occurs for the first few eigenvalues around E — 0. 

Some care must be taken in establishing exactly what gap is to be considered 
in the context of finite-size scaling, since the gap is a function of k — kc, 
that is, how far above or below the ESQPT the gap is measured. The phase 
transition does not fall exactly "on" an eigenvalue (k^ is in general noninteger), 
the gap is varying singularly with k at k^, and the quantum corrections are 
fluctuating most strongly for the first few eigenvalues in the vicinity of kc [33]. 
Therefore, in considering the finite-size scaling at the mean field level, it is only 
meaningful to examine the gap some sufficient number of eigenvalues above 
or below the ESQPT, but nonetheless close enough {\k — kd « N) that the 



' As noted above, different normalization conventions may be encountered for the 
model Hamiltonians. Overall multiplcation of the Hamiltonian by a factor gives 
rise to a superficial difference of unity in the finite-size scaling exponents. 
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Fig. 9. Quantitative comparison of quantum and semiclassical results for the gap, 
including finite-size scaling properties, in the vicinity of the ESQPT (i?~0). Cal- 
culations are for the U(3) vibron model I = states with ^ = 0.5. (a,b) Eigenvalue 
spectrum and its first difference, i.e., the gap, shown as functions of {k — kc)/N 
for N = 100 (open circles) and 1000 (solid circles). The semiclassical result (3.10) 
or (3.11) in terms of the W function (with a = 2.49) is shown for comparison (solid 
curve), (c) Scaling of the gap with respect to N, evaluated at fixed quantum number 
k — kc = 5 relative to the ESQPT, for ^ = 0.3, 0.4, . . ., 0.9. The semiclassical results 
for the scaling (with a = 1.24, 1.92, 2.24, and 2.35) are shown for comparison (solid 
curve). The results of the asymptotic logarithmic expression (3.12), evaluated at 
N = 10^, are also indicated (open triangles). 



scaling appropriate to the ESQPT dominates over the usual A ~ A?" ^ scaling. 

The gap for k — kf. = 5 is plotted as a function of N, for 10^ N < 10^, in 
Fig. 9(c). (The quantity plotted is essentially the gap between the fifth and 
sixth eigenvalues above E = 0, but interpolation is necessary, since k — kc is 
discrete and noninteger in the actual spectra.) Note foremost that the gaps 
for ^ = 0.3 and ^ = 0.9, or for ^ = 0.4 and ^ = 0.8, or for ^ = 0.5 and 
^ = 0.7, converge towards each other for large N . Since S(,^) is symmetric about 
^ = 0.6 [see (3.9)], this demonstrates that the asymptotic behavior depends 
on ^ through S(^), as expected if the properties of the ESQPT are dominated 
by the fuj value (3.8) of the parabohc top of the barrier. The semiclassical 
estimate (3.11) is shown for comparison (using only one fixed a value for each 
symmetric pair of ^ values, for simplicity) and appears to reasonably reproduce 
the finite-size scaling. The results of the simple logarithmic approximation 
A^A K, 27rE{^y/'^/N from (3.12), evaluated at iV = 10*^, are also shown for 
reference. 
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4-3 Order parameters 



Let us now consider the singularity in the order parameter (Nb) (or {N2)), 
which plays a defining role for the ground state QPT. The evolution of the 
order parameter {Nij)k is shown as a function of ^ in Fig. 8(c), again for the 
U(3) vibron model, for the same level {k/N = 0.2) considered in Fig. 8(a). 
This quantity is closely related to the energy plotted in Fig. 8(a), since 



by the Feynman-Hellmann theorem. 



N 



(4.1) 



It is seen that {Nb)k undergoes a dip towards zero at C,=C,c^, which becomes 
sharper and deeper with increasing A^. At the scmiclassical level, one of the 
essential characteristics of the ESQPT was localization of the wave function at 
x — 0, together with vanishing classical velocity (hence, — 0). In coordinate 
form, Nfj ocp^/N^ + [with the coordinate definitions used in (3.2)], so the 
natural extension to the fully quantum description is localization of probability 
with respect to occupation number at ^ 0. The order parameter is shown as 
a function of energy in Fig. 8(d), for the same fixed ^ value (^ = 0.5) considered 
in Fig. 8(b). The "evolution" of properties with respect to excitation energy is 
of necessity discrete, since for finite the eigenvalue spectrum is itself discrete 
[Fig. 8(d) inset]. It is apparent from Fig. 8(c,d) that, while {Nb)k drops towards 
zero at the ESQPT, and the dip becomes sharper and deeper with increasing 
N, {Nb)k is far from actually reaching zero at the finite A^ being considered. 



5 Quantum phases 



So far we have considered the excited state quantum phase transition as a 
singularity in the evolution of the excited state properties rather than as a 
boundary between phases. A central question which arises in connection with 
the ESQPT phenomenology concerns the meaning of "phases" for excited 
states, namely, whether or not the excited states on each side of the phase 
transition can meaningfully be considered to belong to qualitatively distinct 
phases. Of course, in thermodynamics, it is well known that phase transitions, 
in the sense of singularities, do not necessarily imply the existence of distin- 
guishable phases, the liquid- vapor transition in the vicinity of a critical point 
being a classic counterexample. Here we approach identification of phases both 
through indirect measures of the structural properties of the states on either 
side of the ESQPT {e.g., order parameters and spectroscopic signatures) and 
directly through inspection of the wave functions. 
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For the ground state, the "phase" is simply indicated by the value of the 
order parameter {Nii)q (or (A^2)o)- In the large limit, the value of {Ni,)q is 
qualitatively different on either side of .^c, namely, vanishing for ^ < and 
nonzero (growing towards for ^ > In contrast, for the excited states, 

{Ni,)k does not show such a qualitative difference between the two sides of the 
ESQPT. Rather, {Nb)k ^ as the level k crosses the ESQPT but is nonzero 
on either side (Sec. 4.3). Therefore, the expectation value {Nb)k by itself does 
not distinguish two "phases" for the excited states. 

The reason is fundamentally related to the classical limit of the problem 
(Sec. 3). Recall that Nb oc p'^/N'^ + x'^. For the classical ground state, the 
kinetic energy vanishes, and the static equilibrium value for x is simply de- 
termined by the location of the minimum in the potential (2). For excited 
states, such a static quantity no longer provides a suitable measure of the 
phase at the classical level, since excited states (with nonzero kinetic energy) 
are not described by a single equilibrium position. Instead, one must con- 
sider a dynamical definition of phase, taking into account the topology of the 
classical orbits in the phase space [8, 9]. The classical analogue of the "expec- 
tation value" of an observable is its time average over the classical motion, 
{f{x,p)) = r^^ f f{x,p)v{E, x)^^ dx. This is also the semiclassical average with 
respect to the first-order WKB probability density P{x) oc v{E,xy ^, so the 
time average carries over naturally to the quantum expectation value. At the 
quantum level, the consequence of the breakdown of the static definition is 
that the expectation value {Nb)k does not provide an unambiguous measure 
of the phase of an excited state. 

Nonetheless, there arc qualitative changes in the spectrum across the ESQPT. 
In particular, the degeneracy patterns with respect to the angular momentum 
(or seniority) quantum number change from those characteristic of the U(ni)(8) 
U(n2) dynamical symmetry above the critical energy to those characteristic of 
the SO(ni +^2) dynamical symmetry below the critical energy. At the critical 
energy, a rapid rearrangement of degeneracies occurs. This is clearly visible for 
all the models in Fig. 6. The evolution of the eigenvalues for the U(3) model is 
shown in detail in Fig. 10, for more angular momentum values {I < 5) and for 
a larger particle number (iV = 25) than in Fig. 6(b). Note that the like-parity 
states (/ odd or / even) form approximate degenerate multiplets [U(2)-like] 
above E = {) for all while multiplets are composed of all / values [S0(3)-like] 
below this energy. For the Lipkin model, the transition between degeneracy 
patterns is understood from the geometric Hamiltonian, as noted in Ref. [4], 
in terms of degenerate parity doublets below the barrier and lifting of this 
degeneracy above the barrier. For the higher-dimensional models, the change 
in degeneracies at the ESQPT is indicative of the breakdown of the adiabatic 
seperation of rotational and radial vibration degrees of freedom at the critical 
energy. 
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Fig. 10. Correlation diagram for the U(3) vibron model (N = 25), with < Z < 5, 
showing the change in angular momentum degeneracies across the ESQPT. 
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Fig. 11. Probability distributions for the entire spectrum of eigenstates, decomposed 
with respect to the Nf, quantum number, i.e., in the U(2) basis, for the U(3) vibron 
model with I = and AT = 1000. The probability distributions are shown for (a) ^ = 0, 
(b) ^ = 0.2, (c) ^ = 0.5, and (d) ^ = 1. Also shown arc the wave functions for 
individual representative states: (e) above the ESQPT = 0.5, k/N = 0.3), (f) at 
the ESQPT = 0.5, k/N = 0.2), (g) below the ESQPT = 0.5, k/N = 0.1), and 
(h) for the S0(3) dynamical symmetry = 1, k/N = 0.1). 
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For the ground state QPT, the persistence of the degeneracies associated with 
the symmetry hmits as the QPT is approached, in spite of strong symmetry- 
breaking interactions, has been explained in terms of quasidynamical symme- 
try [53]. The quahtative distinction between the states on either side of the 
QPT hes in their forming approximate embedded representations of either 
the U(ni) (8'U(n2) or S0(ni + n2) algebras. (In particular, the phases obtained 
on either side of the ground state QPT have been characterized for the IBM 
in Ref. [54].) We therefore note that it is of considerable interest to deter- 
mine whether or not there is a similar sharp distinction between the states, as 
forming approximate embedded representations of one or the other of these 
algebras, across the ESQPT. 

To consider the question of phases further, let us inspect the structure of the 
wave functions for the U(3) vibron model excited states, as decomposed in the 
U(2) (good Nj,) dynamical symmetry basis. Each density plots in Fig. ll(a-d) 
concisely sumarizes the decomposition the entire spectrum of / = eigenstates, 
for a given value of ^. [A horizontal slice across the plot gives the "wave 
function" of one excited state or, more precisely, the squared amplitudes in its 
decomposition with respect to the U(2) basis. The ground state is represented 
by the bottommost shce.] 

To provide context, first consider the structure of the states when no ESQPT 
is present. For ^ = 0. the Hamiltonian is diagonal in the U(2) basis [Fig. 11(a)]. 
As ^ increases towards 0.2, the ground state critical value, there is a spreading 
of the probability distribution over many neighboring basis states, essentially 
confined to a teardrop shaped region of the Nh-k plot [Fig. 11(b)]. At the 
other limit, = 1, where the S0(3) dynamical symmetry occurs, the proba- 
bility decomposition also follows a regular pattern [Fig. 11(d)]. Here the wave 
functions for the eigenstates with respect to the U(2) basis are known analyti- 
cally [30, 55]. The probability distribution for each state is reflection symmetric 
about Nfj/N — 0.5, peaked at two symmetric extreme values. 

The relevant plot for consideration of the ESQPT is now Fig. 11(c), where 
the probability decompositions of the eigenstates are shown for ^ = 0.5. Indi- 
vidual wave functions are shown in Fig. ll(e-g), with an S0(3) wave function 
[Fig. 11(h)] given for comparison. Below the critical energy (k/N pa 0.2), the 
probability decompositions, taken in aggregate, bear a marked resemblance 
to those obtained for the S0(3) dynamical symmetry in Fig. 11(d). However, 
they are scaled towards smaller iV^ and, in particular, are approximately reflec- 
tion symmetric about a reduced value of Nt, {N^/N !v 0.4). In the immediate 
vicinity of the critical energy, the probability is moderately localized at low 
Nb. Above the critical energy, the probability distribution for each eigenstate 
is again strongly peaked at two extreme values of Nb, but these values are not 
symmetric about a flxed Nh as they are below the critical energy. Rather, their 
midpoint increases approximately linearly with excitation quantum number. 
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as reflected in the linear behavior of {Nj,) above £' = in Fig. 8(b). 

The qualitative distinction between the wave functions below and above the 
critical energy is therefore clearly apparent when the states are viewed in 
aggregate, as in Fig. 11(c). For any given interaction parameter value ^, the 
qualitative distinction also apparently involves reflection symmetry (or lack 
thereof) about some fixed Ni^/N < 0.5. However, the appropriate means of 
constructing a simple measure which allows the immediate characterization of 
the "phase" of a state taken in isolation is not obvious and requires further 
consideration. 



6 Conclusions and outlook 

In the present work, we have seen that the ESQPT phenomena are universal 
to a broad family of two-level models with pairing interactions, including not 
only the s-b models {e.g., Lipkin model, vibron model, and IBM) but also 
the generic two-level bosonic and fermionic pairing models. The properties of 
the eigenvalue spectra (including the quantum gap or level density) and order 
parameters in the vicinity of the ESQPT have been investigated, both at the 
semiclassical level and numerically for the full quantum problem. The finite- 
size scaling properties of the gap have been considered at the mean-field level. 
Qualitative differences between distinct "phases" on either side of the ESQPT 
have also been noted. In the process, both parallels with and differences from 
the conventional ground state QPT have been identified. 

An important aspect is the experimental evidence for ESQPTs. This requires 
the identification of physical systems described by algebraic Hamiltonians with 
^ > for which states with quantum numbers k » 1 can be observed. The 
most promising situations thus far are in molecules described by s-b boson 
models [30] . Further examples are needed to fully understand the experimental 
imphcations of the ESQPT. 

At the theoretical level, several interesting questions remain even for the basic 
two-level models considered so far. Here the main intent was to note the as- 
pects of the ESQPT common to the full family of two-level models. A detailed 
investigation of the specific properties of the bosonic and fermionic two-level 
models with nontrivial degeneracies for both levels is in order [27]. In particu- 
lar, an investigation of the U(ni -|- n2)/[U(ni) ® U(n2)] geometry of the models 
is necessary. While the coherent state analysis has been pursued extensively 
for the s-b models [2, 16, 38, 39, 56], only preliminary use of coherent states has 
so far been made for the generic pairing models [57]. 

Also, the role of integrabihty [58] in the ESQPT must be explored. The SO(n)- 
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invariant Hamiltonian (2.2), or more generally the SO(ni) S0(n2)-invariant 
or Sp(ni) (g) Sp(n2)-invariant Hamiltonian (2.7), is integrable and, moreover, 
leads to a separable and effectively one-dimensional problem (3.1) in the clas- 
sical limit. Application of the quantization condition (3.3), which played a 
crucial role in the semiclassical analysis, is limited to integrable (or approx- 
imately integrable) systems. However, interactions beyond the pure pairing 
interaction are necessary for realistic applications. These interactions destroy 
integrability. They also in general give rise to first-order ground state QPTs, 
for which even the ground state scaling properties have only been partially 
characterized [50, 59]. Possible manifestation of an ESQPT in the spectrum of 
a nonintegrable Hamiltonian is discussed in Ref. [60], but a general theoretical 
foundation for ESQPTs in nonintegrable systems has yet to be developed. 

It would be valuable to bridge the gap between the ground state QPT, where 
the quantum properties scale according to power laws, and the ESQPT, where 
the singularity is logarithmic in nature. Scmiclassically, the connection be- 
tween the two is nontrivial, since the ground state QPT involves a pure quartic 
potential {i.e., no barrier) [Fig. 4(b)], while the ESQPT analysis [Fig. 4(c)] 
requires the classical turning point at the barrier to be well-separated from 
the classical turning point at the outer wall of the well [i.e., a high barrier or, 
conversely, small h'^/{2m) ~ l/iV], so that many states lie below the ESQPT. 
In actual spectroscopic applications, such as to nuclei, often only relatively 
few low-lying states can be observed experimentally. Therefore, the finitc-iV 
precursors of the ESQPT in the intermediate regime, where the ESQPT is 
low-lying in the spectrum, are of special interest. 

The main physical interest, however, lies in possible broad relevance of the 
ESQPT phenomena to various mcsoscopic systems, at least those dominated 
by pairing interactions. In this regard, the analysis must be extended to more 
realistic multi-level pairing models. For instance, the multi-level pairing model 
with equally-spaced levels is of special interest for application to the spectra 
of superconducting grains [61]. Multi-level pairing models can also provide 
a foundation for realistic calculations with the nuclear shell model [62]. The 
models considered in the present work may be constructed as the "infinitely- 
coordinated" hmit [49] of the Ising-type spin- lattice models, i.e., the hmit in 
which all sites interact equally with all others by a long-range interaction. 
It would thus be of interest to examine under what conditions an ESQPT 
may occur in such models for finite-range interactions. (It has recently been 
shown [63] that the fermionic two-level pairing Hamiltonian is related to a 
Bose-Hubbard Hamiltonian [64] by an exact boson mapping, suggesting pos- 
sible apphcation of the ESQPT concept to correlated electron systems or to 
ultracold atoms trapped in optical lattices.) An alternate avenue for exten- 
sion to realistic systems is through coupling of multiple two-level systems as 
subsystems, e.g., the Dicke model [65] for quantum optical systems may be 
obtained as two coupled Lipkin models. 
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A Quasispin and multipole Hamiltonians for the two-level boson 
model 



The two-level boson model [Fig. 2(a)] is characterized by two overlaid alge- 
braic structures: a U(n + 1) structure (n = 2L + 1) arising from the bilinears 
in the creation and annihilation operators and an SU(1, 1) structure involv- 
ing the pairing quasispin operators. The twin algebraic structures provide 
a simple relationship between the multipole Hamiltonian (2.2) and pairing 
Hamiltonian (2.3) for the system. In this appendix, we summarize the rel- 
evant algebraic properties and deduce the explicit relationship between the 
pairing and multipole Hamiltonians, for arbitrary L and for the both possible 
phase choices. The relationship noted for the IBM in Ref. [23] is recovered as 
a special case. 

The well-known bosonic SU(1, 1) quasispin algebra [25] or fermionic SU(2) 
quasispin algebra [26] generators are given, in the convention we adopt here, 

by 



m m m 



and obey commutation relations 

[Sj+,SjJ\ — ^2Sjz [Sjz,Sj+] — +Sj+ [Sjz,SjJ\ — —Sj^, (A. 2) 

where the upper and lower signs apply to the bosonic and fermionic cases, 
respectively. The quasispin z projection Sjz is simply related to the occupancy 
Nj = J2m c]jnCjm, by Sjz — ^{Nj ± Qj), where Qj is the half-degeneracy of level 
j. The squared quasispin 

= ^% ^ K^j+^J- + Sj-Sj+) — Sjz{Sjz — 1) T Sj+Sj^, (A. 3) 

with eigenvalues (S^) = Sj{Sj T l)i is conserved separately for each level by 
the Hamiltonian (2.3). Eigenstates are therefore characterized by seniority 
quantum numbers Vj = 0, 1, . . ., defined by Sj = ^{fij ± Vj). 
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For the two-level bosonic system defined in terms of a singlet level s^'^^ and an 
n-fold degenerate level b^^^ (n = 2L + 1), the quasispin generators (A.l) are 
given explicitly, in tensor notation, by 



= + =i(iV, + L + i), 



(A.4) 



where = (-)^-'^Tl^j and C/W • T/W = (-)^(2L + iy/^{A x = 

E^(-)^f/i^^1^iJ . The generators of the SU,(1, 1) and SU6(1, 1) algebras can be 
combined to form a sum-quasispin algebra with two possible relative phases, 
yielding SU±(1, 1) algebras with generators 

SU±(1,1): S^ = Ss+±S,+ S. = Ss-±S,. S,^Ss, + Sbz. (A.5) 

The subalgebras of SUs(l, 1) (8) SUb(l, 1), and their associated quantum num- 
bers, are 

SU+(1,1) 



v+ 



SU,(1, 1) SUfe(l, 1)D { 

Vs=0, 1 "6 



Us(l) ® U5(l) (A.6) 
SU_(1,1). 



The SUs(l, 1) algebra for the singlet s-boson level is trivial, in that = —3/4 
identically by application of the canonical commutation relations for Sq and 
So- This constrains Vs to the values or 1. Since Sz — S is integral for SU(1, 1) 
representations, it follows that Vs = for A'^s even and Vs = l for Ns odd {i.e., 
Vs = Ns mod 2). 

At fixed total particle number N (and therefore fixed Sz), the operator S'+S'_ 
is trivially related to the Casimir invariant of SU±(1, 1) by (A. 3), as 
S+S- = S^{Sz - 1) - S^, with eigenvalues \[N{N + 2L) - v±{v± + 2L)]. A 
pairing Hamiltonian (2.3) with pairing interaction chosen proportional to the 
SU±(1, 1) Casimir operator is thus given by 

{Hpp)± = sN, + 4«:(-)^+i(^,+ ± S,+){Ss- ± S,_), (A.7) 

where the coefficient on the last term is chosen for convenience below. 

The two-level boson system is alternatively characterized by the Lie algebra 
U(n -|- 1), with tensor-coupled generators 

U(n+ 1) : {s^ X §p {s^ x (6^ x 5)^^) {b^ x b)^^\ (A.8) 
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for A = 0, 1, ... , 2L. Two distinct SO±(n + 1) subalgebras are obtained, with 
generators 



SO±(n + 1) : 



i[(st X - (6t X 



(fet X (A.9) 



for A restricted to odd values. In the case of the IBM, it is the S0+(6) alge- 
bra which contains the physical quadrupolc operator [66]. [Therefore, conven- 
tionally, the S0_|_(6) algebra is simply denoted by S0(6), while the alternate 
S0_(6) algebra is denoted by S0(6).] An n-dimensional rotation algebra SO(n) 
is obtained by retaining only the generators {h'^ x h)^^^ (A odd), and an S0(3) 
algebra by retaining only (6^ x 6)^^^ [19]. 

The subalgebras of U(n -|- 1), and their associated quantum numbers, are thus 

' SO+(n + l) 



U(n+ 1)D < 

N 



0-+ 

U,(l) ® U(n) 

SO_(n + 1) 



> DS0(n)DS0(3) 

V J 



(A.IO) 



The Casimir operators of the subalgebras are, explicitly, 

C2[S0±(n + 1)] = (±)2(st xb±b^ X sY^^ ■ (s^ x 6 ± fet x + C2[S0(n)] 
C2[S0(n)] = 4 ^ (fet X 6)W ■ (fet X 6)^ 

^ ~ (A 11) 

C,[V{n)] = i-)'^ib^-b) = N, ^ ■ ^ 

C2[U(n)] = ^(6^ X • (6t x 



with eigenvalues 2a±{a±+2L), 2v{v+2L—l), Nf, and N{,{Nh+2L), respectively. 
The mathematically natural Casimir form of the Hamiltonian is 



K 



iHc)±^eN, + -C2[S0±in + l)]. 



(A.12) 



The two phase choices for SO-(-(n + 1) yield identical eigenvalue spectra but 
different eigenstates [66]. For physical reasons, the SO(n)-invariant "multipole- 
multipole" Hamiltonian 



MM)± 



eNb + /«(±)[(s^ xb±b^ X s)^^^ ■ {s^ xb±b^ X s)^^^] (A. 13) 



is commonly used. At fixed v, He and Hmm differ only by a constant offset, 
with {Hmm)±^{Hc)± - («:/2)C2[SO(n)], by (A.ll). 

To relate the U(n -|- 1) and SU(1, 1) descriptions, let us first observe that the 
Casimir operators of SO(n) and SU(,(1, 1) are related, and that the SO(n) an- 
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gular momentum quantum number v and the SU6(1, 1) seniority Vb are actually 
identical. This is an example of a general correspondence between the algebras 
SO(n) and SU(1, 1) [67,68]. According to the basic quasispin relations above, 

4S^ = (TVft + L + i)(7V6 + L - I) - ASb+Sb-, (A. 14) 

and 4:Sb+Sb- = (fe^ • b^)(b ■ b) can be recast in terms of coupled bilinears in the 
creation and annihilation operators {e.g., by the tensor contraction methods 
of Ref. [69]) as 

4Sb+Sb- = (-)^+'(&^ • b) + ^{-)^{b^ X 6)(^) • {b^ X b)^^\ (A.15) 

A 

The right hand side we reexpress in terms of the Casimir operators (A. 11) as 

4Sb+Sb- = -Nb + C2[U(n)] - |C2[S0(n)]. (A.16) 

On the other hand, (4S^) = {vb + 2L + 1) (f + 2L — |) from the definition of se- 
niority. Comparison with the result for (4S^) obtained from (A. 14) and (A.16) 
establishes the identity of v and Vb- 

Comparison of the full Hamiltonians Hpp (A. 7) and Hmm (A. 13) is then 
straightforward, with the aid of (A.16), yielding 

{Hmm)± = {Hpp)i^i-)L] + k{-)''[N{N + 2L) - v{v + 2L- 1)]. (A.17) 

Thus, for L even, the Hamiltonians {Hmm)± and (Hpp)^ differ only by a 
constant {i.e., a function of conserved quantum numbers), while, for L odd, it 
is {Hmm)± and {Hpp)± which differ only by a constant. Observe, therefore, 
that for the U(2), U(6), etc., models {L even), the conventional SU+(1, 1) phase 
choice for the pairing interaction {Gij — 1) corresponds to the unconventional 
SO_(n + 1), or SO(n + 1), phase choice for the multipole interaction, and vice 
versa. 

Finally, let us consider how the general U(n + 1) scheme (A. 10) specializes 
to the Schwinger realization of the Lipkin model {L = 0). With the usual 
Schwinger angular momentum operator definitions J+ = b^s, J_ = s^b, and 
Jz = ^{b^b — s"^s), the subalgebra structure is 



U(2)dSU(2)d < 



S0+(2) = S0,(2) 



U(1) = S0.(2) (A.18) 
S0_(2) = S0y(2) . 

Jy 

The breaking of U(n) into S0+(n + 1), S0_(n + 1), and U(n) subalgebras 
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therefore has a particularly simple interpretation, as the projection of the 
SU(2) angular momentum alternatively along the y, or z axes. 



Furthermore, both the s and b bosonic levels are singlet levels in the Lipkin 
model. It follows, as already noted for the s boson in (A. 6), that the seniorities 
Vg and Vf, are restricted to the values or 1, with Vg = Ng mod 2 and Vf, = Nh 
mod 2. Thus, the grading quantum number g defined in Sec. 4.1, which deter- 
mines the Lipkin model parity n = (— )^, is simply the SU6(1, 1) quasispin se- 
niority Vb. This formally explains the similarity, discussed in Sec. 4.1, between 
the role of g in the Lipkin model and that of the SO(n) angular momentum 
in the higher-dimensional s-b algebras [Fig. 6(a,b)]. However, for the Lipkin 
model (n = 1), there is no SO(n) angular momentum dual to the SUb(l, 1) 
seniority. 



B The Lambert W function 

The Lambert W function [70] is implicitly defined as the solution y = W{x) 
to the equation 



In this appendix, we summarize the essential properties of the W function 
needed for the present analysis (Sec. 3). The complex analysis, asymptotics, 
series expansion, etc., of the W function are considered in detail in Ref. [70]. 

Considered as a real- valued function of a real variable, W{x) is single- valued 

for a; > but double- valued for — 1/e < ,x < 0, with branches Wo{x) > —1 and 
W^i{x) < — 1. These branches are plotted in Fig. B.l. The function has the 
asymptotic form [70, 71] 



as a; — >• , shown as the dashed curve in Fig. B.l. 

From the defining equation (B.l), it follows that W obeys the identity 



X — ye^. 



(B.l) 



W^i{x) ~ log(-x) - log[-log(-x)] 



(B.2) 



\ogW{x) 



logx — W{x). 



(B.3) 



Differentiation yields 



W'{x) 



W{x) 



(B.4) 



It also follows from (B.l) that the 



x[i + w{x)y 

equation 



y\ogy + cy ^x 



(B.5) 



has solution y = x/W{e'^x), as needed for Sec. 3.3. 
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X 

Fig. B.l. The Lambert W function. For — l/e<x<0, W{x) is double- valued, with 
branches Wq{x) > — 1 and W-i{x) < —1. The asymptotic form of W-i{x) given 
by (B.2) is shown for comparison (dashed line). 
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